Interplay between interaction and chiral anomaly: Anisotropic perfect metallicity in 
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We uncover that long-range interactions give rise to perfect metallicity (zero resistivity) at zero 
temperature in Weyl metals when electric fields are applied to the direction ol the momentum 
vector to connect two paired Weyl points. Performing the renormalization group analysis, we find 
that the distance between two Weyl points becomes enhanced logarithmically at low temperatures 
although the coupling constant of such interactions vanishes inverse-logarithmically. Considering 
the Adler-Bell-Jackiw anomaly, scattering between these two Weyl points becomes suppressed to 
increase electrical conductivity in the "longitudinal" direction, counter-intuitive in the respect that 
interactions are expected to reduce metallicity. Correlations with topological constraints allow 
unexpected and exotic transport properties. 
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Nontrivial global structures of ground states some- 
times violate classically respected conservation laws at 
quantum levels, referred to as (quantum) anomalies [1]. 
When anomalies are associated with breaking of local 
(gauge) symmetries, it means that their corresponding 
quantum theories are not consistent and such anoma- 
lies should be cancelled, introducing meaningful quantum 
fields. Indeed, the standard model and string theories 
are constructed consistently, cancelling gauge and gravi- 
tational (also conformal) anomalies, respectively [2]. On 
the other hand, when such anomalies are related with 
breaking of global symmetries, they give rise to vari- 
ous interesting physical properties. In particular, vari- 
ous types of topological terms associated with quantum 
anomalies arise to play essential roles in quantum criti- 
cality of quantum matter [3] . In addition, they turn out 
to be responsible for quantum number fractionalization, 
given by Goldstone-Wilczek currents [4]. Actually, an 
emergent non-abelian chiral anomaly has been proposed 
to cause so called deconfined quantum criticality in low 
dimensional spin systems [5, 6]. Furthermore, such topo- 
logical terms sometimes give rise to anomalous (quan- 
tized) electrical or thermal Hall effects [7, 8]. Quantum 
anomalies govern quantum criticality, quantum number 
fractionalization, and anomalous transport phenomena 
[9]. 

In this study we focus on the role of the Adler- 
Bell-Jackiw anomaly or chiral anomaly [1] in anoma- 
lous transport phenomena. This anomaly means that 
classically conserved chiral currents, that is, currents of 
right-handed (Weyl) fermions minus those of left-handed 
fermions, are not preserved at quantum levels due to 
nontrivial global configurations of gauge (or electric and 
magnetic) fields. Applying magnetic fields to gapless 
semi-conductors, a Dirac point described by the four- 
component Dirac spinor splits into two Weyl points gov- 
erned by the two-component Weyl spinors with opposite 
chiralities, where the distance between two Weyl points is 



proportional to the applied magnetic field [10] . The chiral 
anomaly gives rise to a topological constraint in dynam- 
ics of Weyl fermions, where right-handed Weyl fermions 
at one Weyl point should scatter into left-handed Weyl 
fermions at the other Weyl point, when currents are 
driven to the same direction as the momentum to connect 
these two Weyl points [11]. Even if short-range scatter- 
ers are taken into account, scattering between these two 
Weyl points becomes suppressed due to the finite dis- 
tance in the momentum space. As a result, the longitu- 
dinal (E || B) magneto-conductivity is enhanced, which 
turns out to be proportional to the square of the applied 
magnetic field or the distance of two Weyl points [11, 12]. 

In this letter, we investigate effects of interactions on 
the "longitudinal" "magneto" -transport in Weyl metals. 
Here, " " will be clarified later. It is almost trivial to ob- 
serve that local four-fermion interactions are irrelevant 
at low energies in a perturbative sense since the density 
of states vanishes at zero energy. Long-range Coulomb 
interactions have been investigated both extensively and 
intensively for transport phenomena in graphene [13]. In 
addition, transverse gauge interactions have been also 
discussed in Weyl- or Dirac- type systems [14]. Recently, 
the chiral anomaly has been calculated in the Weyl sys- 
tem [15]. However, the interplay between long-range in- 
teractions and the chiral anomaly has not been investi- 
gated clearly. In particular, it remains mysterious how 
this combination gives rise to anomalous "longitudinal" 
"magneto" -transport phenomena. 

Performing the renormalization group analysis, we re- 
veal that the distance between two Weyl points becomes 
enhanced logarithmically at low temperatures although 
the coupling constant for transverse long-range interac- 
tions vanishes inverse-logarithmically (expected in three 
dimensions). This is in contrast with "conventional" 
Weyl metals without interactions [10-12], where the dis- 
tance between two corresponding Weyl points remains 
finite. As a result, scattering between two Weyl points 
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becomes suppressed much more than the case of nonin- 
tcracting Wcyl metals, which increases electrical conduc- 
tivity in the direction to connect the momentum vector 
between two Weyl points. We argue that anisotropic per- 
fect metallicity arises in the case of a finite chemical po- 
tential, where perfect metallicity appears for the longitu- 
dinal direction while normal metallic behaviors result for 
other directions. We discuss this interaction-enhanced 
metallicity based on the quantum Boltzmann equation 
approach in the presence of both long-range interactions 
and the chiral anomaly. 

We start from quantum electrodynamics with a topo- 
logical 9-term in three spatial dimensions (6— QED 4 ) 

£=-^"" + ^ + 0^'-', (1) 

where tp is a four-component Dirac spinor to take both 
chirality (associated with either orbital or sublattice in- 
dices) and spin quantum numbers and is an electro- 
magnetic vector potential regarded as a quantum field. 
Dfj, = + ieAp in if) = j^D^ is a covariant derivative 
with an electric charge e, where 7^ is the Dirac gamma 
matrix satisfying the Clifford algebra with \i = 0, 1, 2, 3. 
pnv _ L_ e nvpa j s a ma g ne tic dual tensor of the elec- 
tromagnetic field strength tensor. One may consider that 
this field theory results from certain lattice models with 
spin-orbit interactions for topological insulators except 
for long-range transverse interactions [16]. However, we 
would like to emphasize that even gauge fluctuations can 
emerge in some lattice models, supporting so called topo- 
logical spin liquids [17]. The 8— term is the fingerprint of 
the topological insulator in three dimensions, the source 
of the (longitudinal) magnetoelectric effect or equiva- 
lently, the half-quantized Hall conductance on its surface 
[8]. 

The Adler-Bell-Jackiw anomaly states that the classi- 
cally conserved chiral current is not conserved at quan- 
tum levels [1], given by 
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with the chiral current = ^7^75^ — -^f — J fit men- 
tioned before, where the 75 matrix is a Dirac matrix 
that anti-commutcs with other Dirac matrices. In other 
words, when electric fields are applied in parallel with 
magnetic fields, the chiral current is not conserved. It is 
important to realize that the 0—term is a boundary term, 
implying that the coefficient cannot be renormalized by 
interactions. However, if inhomogeneous magnetic fields 
can be applied to this topological insulating state, the 
9 coefficient depends on position [16, 18]. As a result, 
this term is not a boundary term any more, which can 
be renormalized by interactions. 

Resorting to this anomaly equation, we rewrite the ef- 
fective field theory as follows [19] 



C = - l -F lxv F^ 



where c M with the 75 Dirac matrix is a chiral gauge field, 
given by d^O. It is interesting to observe that the Dirac 
point splits into two Weyl points at K = ±c. In this 
respect our problem is to investigate the nature of the 
quantum critical point between a topological insulator 
and a band insulator in the presence of inhomogeneous 
magnetic fields. In other words, we study how both the 
interaction parameter e and the distance between two 
Weyl points are renormalized to evolve at low tem- 
peratures, and reveal how these renormalization effects 
modify transport properties, compared with the case in 
the absence of interactions. Such inhomogeneous mag- 
netic fields may be created by either ferri-magnctism [16] 
or some ferromagnetic clusters, given by randomly dis- 
tributed magnetic ions [20]. 

Introducing counter terms, we rewrite this effective 
field theory as follows [19] 



C = ~F ltv F'"' + ^(Z^ 



- Z c fa 5 - Z e e4) ip, 
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where Za, Z c , and Z e are field renormalization con- 
stants of ip, A^, c M and the vertex or coupling (e) renor- 
malization constant, respectively. We emphasize that the 
renormalization factor Z c has never been introduced as 
far as we know, thus regarded as an essential aspect of 
this study. We alert that c M is a background gauge field, 
not dynamical in the present study. 

In order to perform the renormalization group analy- 
sis, we resort to the dimensional regularization. A subtle 
point arises due to the presence of the 75 matrix, which 
needs some care for its treatment, because its existence 
depends on dimensionality [1]. We point out that the 
presence of the 75 matrix makes our calculations much 
more complicated and laborious. One nontrivial check for 
the validity of our calculations is that the Ward identity 
is respected in the one-loop level. All details are pre- 
sented in our supplementary material [19]. As a result, 
we obtain our coupled renormalization group equations 
for e and c M , 
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The first equation is nothing but the conventional 
renormalization group equation of the coupling constant, 
which tells us that the electric charge renormalizes to 
vanish at zero temperature. Even if the coupling con- 
stant vanishes, the background chiral gauge field flows to 
go to infinity. Inserting the solution of the first equation 
into the second equation, we find [19] 
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where cJ^ and e u are chiral gauge field and coupling con- 
stant at the energy scale of the bandwidth or cutoff. The 
chiral gauge field increases in a logarithmic way. This 
indicates that the distance between two Weyl points be- 
comes infinite at zero temperature, implying that their 
scattering events are suppressed completely. This renor- 
malization effect should be observed in transport coeffi- 
cients. 

Our framework for anomalous transport is the "semi- 
classical" quantum Boltzmann equation approach. Here, 
the term "semi-classical" means that the role of both 
Berry curvature and chiral anomaly or the topological 
9— term is introduced from coupled semi-classical equa- 
tions of motion based on the wave-packet picture in solids 
[7]. We incorporate this information into the quantum 
Boltzmann equation, where the case with a finite chem- 
ical potential is considered, more generic than the case 
with two Weyl points. In principle, one can derive the 
quantum Boltzmann equation in a matrix form, regarded 
as a full quantum transport theory [21]. However, its 
derivation is much complicated and not easy to per- 
form. We would like to mention that our phenomcnolog- 
ical "quantum" transport theory is to generalize a recent 
"classical" Boltzmann transport theory [12]. As a result, 
inelastic scattering events can be taken into account nat- 
urally in the presence of the topological 6— term. 

Based on our phenomenological quantum Boltzmann 
equation approach, we propose the following expression 
for the longitudinal conductivity [19] 

a L (T)^(l + /C[c(T)] 2 K(T), (9) 

where K. is a positive numerical constant and the normal 
conductivity a n (T) is determined by (gauge- interaction 
induced) intra-scattering events at one Weyl point. An 
essential point is in the c(T) term, where c(T) is the dis- 
tance between two Weyl points, given by our renormal- 
ization group analysis. Here, the vector index p is fixed 
and omitted for simplicity. Actually, this expression is to 
replace the applied magnetic field in noninteracting Weyl 
metals with the distance between Weyl points in inter- 
acting Weyl metals, when electric fields are applied in 
parallel with "magnetic fields" or the momentum vector 
to connect such Weyl points, implying the reason why 
we call "longitudinal" in front of the conductivity. We 
emphasize that this expression recovers that of the origi- 
nal proposal, where the magneto-conductivity is propor- 
tional to the square of the distance between two Weyl 
points. Although the distance does not renormalize in 
the noninteracting case, it increases logarithmically in 
the presence of long-range interactions. 

One aspect should be pointed out carefully. We proved 
that the distance between two paired Weyl points in- 
creases logarithmically to infinity at zero temperature, 
originating from long-range transverse interactions. It 
should be noticed that this result appears at zero chem- 
ical potential. On the other hand, we found the longi- 
tudinal conductivity [Eq. (9)] for a finite chemical po- 
tential. Can we expect the similar enhancement of the 



distance between two paired Weyl points in the case of 
a finite chemical potential? We would like to point out 
that it becomes more delicate to perform the renormal- 
ization group analysis in the case of a finite chemical 
potential because existence of abundant soft modes near 
Fermi surfaces can give rise to various types of competing 
instabilities. However, we expect that the renormaliza- 
tion group equation [Eq. (6)] for the chiral gauge field 
will hold at least in the one-loop level as long as the re- 
sulting state remains "normal" against symmetry break- 
ing. Then, the positive slope with respect to the scale 
parameter In p in Eq. (6) indicates that the background 
chiral gauge field increases monotonically as temperature 
decreases down to zero. On the other hand, the curva- 
ture in the curve of the chiral gauge field as a function of 
the scale parameter turns out to be proportional to the 
beta function of the coupling constant. As long as the 
beta function of the coupling constant is positive, we con- 
clude that the distance between two paired Weyl points 
increases to infinity at zero temperature. Actually, it is 
straightforward to see that the gauge coupling constant 
decreases monotonically to vanish at zero temperature, 
implying the positiveness of its beta function, where the 
presence of Fermi surfaces weakens long-range interac- 
tions between electrons more drastically than the case in 
the absence of Fermi surfaces [22]. 

We plot the longitudinal resistivity in Fig. 1, given 
by the inverse of the longitudinal conductivity [Eq. (9)]. 
The normal resistivity is 

p n {T)=p +TT a , 

where po is a residual resistivity coming from disorder 
scattering and J 7 is a positive numerical constant. The 
exponent a is 5/3 in the three dimensional fermion-gauge 
coupling problem with a Fermi surface [14]. But, we 
would like to mention that the exponent can be regarded 
as a fitting parameter when there are other types of elec- 
tron correlations. An essential point is the temperature 
dependence of the background chiral gauge field. Since it 
diverges logarithmically, we reach a remarkable conclu- 
sion that the resistivity vanishes identically at zero tem- 
perature, regardless of disorder scattering. This is rather 
unexpected because interactions give rise to localization 
conventionally instead of derealization. The interplay 
between the chiral anomaly and long-range electron cor- 
relations can cause anomalous transport phenomena. 

However, it should be noted that this perfect metal- 
lic behavior arises only when electric fields are applied 
to the direction of the momentum vector to connect two 
paired Weyl points. Resistivities of other directions are 
quite conventional. Strong anisotropy in transport is an 
important character of Weyl metals. Unfortunately, this 
perfect metallic behavior may be difficult to observe be- 
cause there exist additional charge carriers from other 
bands at a finite chemical potential [23]. 

In summary, we investigated the quantum critical 
point of the topological phase transition from a topolog- 
ical insulator to a band insulator in the presence of in- 
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FIG. 1: Temperature dependence of longitudinal resistiv- 
ity in an interacting Weyl metal. We plot a dimension- 
less longitudinal resistivity scaled by a residual resistivity 

po, = (l + ZCciJl + lllnfiy'll + ^/po)^] as 

a function of temperature T of a unit of K, using D = W 4 K, 
K. = 0.1, cd — it, T/po — §.\K~ a , and a = 5/3. In particu- 
lar, we used the value of eu that reproduces the fine structure 

constant, i.e., -At = 1/137. An essential feature is the emer- 
gence of perfect metallicity in the direction of the momentum 
vector between two paired Weyl points. See the inset figure, 
which plots the same figure but in a low temperature region, 
where the thick (red) line denotes the longitudinal resistivity 
of an interacting Weyl metal and the dashed (green) line rep- 
resents that ( EI jp- = 1 + ^■T a ) of a "normal" metal with 
gauge interactions. Generally speaking, metallicity becomes 
enhanced in Weyl metals as a result of the interplay between 
transverse long-range interactions and the chiral anomaly. 



homogeneous ferromagnetism or under nonuniform mag- 



netic fields, where the topological 6— term gives rise to 
a topological constraint in dynamics of bulk fermions, 
referred to as the Adler-Bell-Jackiw anomaly. Such inho- 
mogeneous magnetic fields serve background chiral gauge 
fields, splitting the Dirac point into two Weyl points. 
Introducing long-range transverse interactions and per- 
forming the renormalization group analysis, we uncov- 
ered that the distance between these two Weyl points 
becomes enhanced to flow to infinity at low temperatures 
although the coupling constant vanishes as expected. Re- 
sorting to the semi-classical quantum Boltzmann equa- 
tion approach, we claimed that the enhancement of the 
distance strengthens metallic properties at low tempera- 
tures when electric fields are applied to the same direction 
as the momentum to connect these Weyl points because 
scattering between the Weyl points are suppressed due 
to their huge distance in the momentum space. 

There remain two important problems. The first ques- 
tion is what happens if we take into account chiral gauge 
fields quantum mechanically. This situation arises when 
ferromagnetic phase transitions occur near the topologi- 
cal phase transition [16]. Is it possible to obtain a novel 
interacting fixed point? The second question is more 
practical thus experimentally verified. If we introduce 
weak anti-localization corrections in the transport theory, 
how is the longitudinal transport coefficient modified? 
This question is still meaningful even without interac- 
tions because this transport signature can be measured 
actually [23]. 

KS was supported by the National Research Founda- 
tion of Korea (NRF) grant funded by the Korea govern- 
ment (MEST) (No. 2012000550). 



Appendix A: Renormalization group analysis in the presence of chiral anomaly 



1. Introduction of counter-terms 



We start from the following Lagrangian 

C B = i^B'fd^B - egAg^B^B - -Fg^Fg" + ^^ e ^Fg^Fg pa , (Al) 

where ipg, Ag^, eg, 9g are bare quantities to be renormalized. Introducing renormalization factors of tpg — 

1/2 1/2 1/2 

Ag p = Afi, egZ A Z^ — Z e e, 9g = Z c 8, one can rewrite the above bare Lagrangian in terms of its 
renormalized part and counter-term part, 

£-b = C r + £ c .t. 

C r = i^d^ - eA^i> - \f^F^ + ^L e ^F^F pa 

Cct. = S^rd^ - 5 e eA^ri> ~ 5 -fF pu F^ + 5 c -^e^ prT F^F pa , (A2) 

4 1D7T 

where 6$ = Z^ - 1, S e = Z e - 1, S A = Z A - 1, 6 C = Z c - 1. 

Incorporating the anomaly equation [Eq. (2)] for renormalized fields into the above expression, we obtain Eq. (4) 
for our renormalization group analysis. This procedure can be performed in a more formal way. This Lagrangian 
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functional is invariant under the chiral transformation of tp — > e ta ^ l5 ip as long as fermions remain massless. Consider 
the following replacement 

Z = j DAD^D^ exp ji J d 4 xC r + C c . t . j 

-> Z = J DAD^D^ exp ji ^ d 4 x£ r + £ c . t . + a(x) J 75 " + i^ 6 ^^^) } . ( A3 ) 

where gauge fixing is assumed. a(x) is an arbitrary and infinitesimal local parameter. Taking a(x) = —9(x), we 
see that the OFF term is replaced with the chiral gauge- field term in Eq. (4). As a result, we obtain the following 
expression 

Z = f DADtpDijj exp j, f d^xC, + 

A.c.t. = ~F^F^ + i^y"^ - ^75)^ - ^eA^V, (A4) 

where gauge fixing is also assumed. It is important to notice that the chiral anomaly equation is satisfied for 
renormalizcd fields, not bare fields. 



2. One-loop structure of QED4 with a background chiral gauge field 

We perform one-loop renormalization group analysis in the presence of the background chiral gauge field, i.e., c M 
= constant. We obtain renormalization constants from £i(f>, c), F^(p,p' , c), U.^(q, c), corresponding to one-loop 
fermion self-energy, one-loop vertex correction, and one-loop gauge-boson self-energy, respectively. 

It is important to notice that the background chiral gauge field is taken into account non-perturbatively. In other 
words, our vacuum state is a state that chiral currents are flowing. Thus, the fermion propagator is modified to be 

1 _ _L_ _ ( P 2 + c 2 )f + 2(p ■ c)f l5 - (p 2 + c 2 )^ 75 - 2(p ■ c)£ 

f p+fij5 (p- C) 2 (p + C) 2 

In this respect the key point is how self-energies of fermions and gauge bosons and vertex corrections are modified in 
this novel vacuum state. A subtle point arises due to the presence of the 75 matrix in the regularization procedure 
[1]. When dimensional regularization is used to regularize loop- integrals including 75, some anomalous terms appear. 
They originate from components perpendicular to physical four-dimensions. We separate out these perpendicular 
momentum components of 1^ = 1^ — W explicitly in our dimensional regularization. However, it turns out that they 
do not result in divergent contributions. 
First, we calculate the fermion self-energy 

-iE,r# c) (- ie ) 2 f ddk (P ± c2)t 2(fc ' ± 2(fc ' c)h5 (fc2 ± c ^ 75 -,-» 

i^,c)-{ ze) J (2 ^ )d7 (fc + c)2(fc _ c)2 7 {p _ k) 2 

= - p2 1 (trip-knk-cnk+c) 2 < (fc2 + c2)rh » 2{k ■ c)7 ^ + 2(fc ' c)7 ^ 757 - {k2 + c2 ) 7 ^ 757 -} ■ 

(A6) 

Replacing the denominator with Feynman parameters, we rewrite the above expression as follows 

f 1 f d d l 1 

-iEiQ*,c) = - 2e2 J o dxdydz 5{x + y + z - 1) J -^-^ — - ^ 3 

x {-2((Z + a) 2 + c 2 ){t + $ + 4(i • c + a ■ c)£ + 4(Z ■ c + a- c)(/ + <f) lb - 2{{l + a) 2 + c 2 )£ lb ) , 

(A7) 

where 1^ = k^ — a p , a p = xp^ + (y — z)c^, A = (x 2 — x)p 2 + (y 2 — y + z 2 — z — 2yz)c 2 + (2xy — 2xz)(p ■ c). Note 
that the degree of divergence in each term depends on only the power of redefined loop momenta Since we need to 
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calculate only divergent terms, we consider 



A] 3 



f 1 f d d l 1 

-iSi(j*,c) = - 2e y o dxdydz 5(x + y + z-l) J ^-y d — - 

x {-2/V - 4^ra A7K + 4Z A Z K c A 7 K 75 - 2/ 2 ^} + /iraie. 

= -2e 2 / dxdydz S(x + y + z — 1) x 
Jo 



d d Z Z 2 \ / /" d d Z l x l K \ 1 

-% + ^5)( / ( 27r)rf [Z 2_ A] 3 j- 4 K7 K -CA7.75)^ (27r)rf [p _ Ap J } + /2mte - 

(A8) 



Integrating over the loop momenta Z, we obtain 



e 2 i f 1 (3 1 1 

-iSi(^, c) = — / dxdydz 8(x + y + z - 1) < -ft + -^75 > + finite 

„2„- 

^75) + finite (A9) 



87r 2 e 
with e = 4 — rf. 

Second, we calculate the vertex correction in the same way as the above 

/■ tf^fc f ^(fc^ + c 2 )//' -2(fc'.c)^ + 2(fc'- c)//' 7 5 - (fc' 2 + c 2 )^75 



(27T) d t (fc'-c) 2 (fc' + C ) 2 



_. 2 f cZ d fc 

" ~ ie J w?' 



„ . (fc 2 + c 2 )^ - 2(fc ■ c)^ + 2(fc ■ c)# 75 - (fc 2 + c 2 )^ 75 p \ -ir, vp 

X7 * (fc- C ) 2 (fc + c) 2 7 J (fc-p) 2 



1 



(fc - p) 2 (k> - c) 2 (fc' + c) 2 (fc - c) 2 (fc + c) 2 

x 7 "{(fc' 2 + c 2 )k" - 2(fc' • c)^ + 2(fc' • c)fe*75 - (fc' 2 + c 2 )^} 

x 7 ^{(fc 2 + c 2 )# - 2(fc • c)£ + 2(fc • c)# 75 - (fc 2 + c 2 )^}lu, (A10) 

where p 2 = p' 2 = m 2 = 0, k^ + q^ = fc^, p^ + q^ = p'^. Similarly, replacing the denominator with Feynman parameters 
and taking only divergent terms, we obtain 

f 1 f d d k 1 

T^(p,p' , c) = — 24ze 2 / dxdydzdudv 5(x + y + z + u + v — 1) 



(27r) d [Z 2 - A] 5 

X7 ,y {((Z + a + g) 2 + c 2 )(/ + ^ + d) - 2(1 ■ c + a ■ c + q ■ c)fi + 2(1 ■ c + a ■ c + q • c)(/ + (fi + ^75 - ((Z + a + 9 ) 2 + c 2 )^} 

x 7 ^{((Z + a) 2 + c 2 )(/ + </l) - 2(1 ■ c + a ■ c)£ + 2(Z • c + a • c)(/ + ^75 - ((Z + a) 2 + c 2 )^^ 

f 1 ( f d d k (l 2 ) 2 l\l 1 

= 48ie 2 J dxdydzdudv 5(x + y + z + u + v - 1) |y ^-^ ^ 2 _ ^ P 5 j {7 A 7 M 7 P } + /mite. 



e 2 



7^ + /imte, (All) 
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where Z M = fc M - o M , = xp p - (y + z)q ll + (y - z + u - u)c M . 

Third, we evaluate the vacuum polarization tensor. Although the procedure is essentially the same as before, this 
calculation is much more complicated, given by 

" dh r M , (fc 2 + c 2 )$ - 2(k ■ c)fi + 2(fc ■ c)j&75 - ( fc2 + c 2 )^7s 
7 1 (k-c) 2 (k + c) 2 

v ((fc + q) 2 + c 2 )(^ + d) - 2(k ■ c + q ■ c)t + 2(fc • c + q ■ c)(ft + ^75 - ((fc + Q? + c 2 )^75 
X7 l- 



l U^(q,c) = -(- l e) 2 J^ [ Tr 



-e 2 



(k + q-c) 2 (k + q + c) 2 

d d k 1 

(2ir) d (fc - c) 2 (fc + c) 2 (fc + q - c) 2 (k + q + c) 2 
x Tr [^{(k 2 + c 2 )^ - 2(fc • c)^ + 2(fc • c)# 75 - (fc 2 + c 2 )^ 75 } 

X7"{((fc + 9) 2 + c 2 )(^ + ^) - 2(fc • c + g • c)^ + 2(fc ■ c + q ■ c)(t + ^)75 - ((fc + ?) 2 + c 2 )^} ] 

(A12) 
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Notice that the extra (—1) factor comes from the fermion loop in the diagram. Now, we have to change the loop 
momenta as before. Straightforward but rather tedious algebras give us the following expression 



iU^ v (q, c) = -6e 2 J dxdydz S(x + y + z - 1) J 



d d l 1 



(2ir) d [I 2 - A] 4 

x[ {i 2 f{-^} + 8{i 2 ) 2 i^i v 

+ (Z 2 ) 2 {4a^ + 4a"6" - V"> • b) - 4ie^ pa c p q a + 8c"c" - \2rf v c 2 - tof (a 2 + b 2 )} 
+l 2 Fl x {8h x h" - 32c x c v } + l 2 n x {8h x h" - 32c x c"} 

+l 2 l x l K {~8h x h K n^ - l6a x b K n» v + 16c x c K rT + mc x b a n Kp e^° + mc x a p n Ka e^} 
+l 2 m u {8a 2 + 8b 2 + 16c 2 } + Fl u l x l K {2,2a x b K + 32c K c A } ] 
+ finite, (A13) 



where 



^ k^ a,^, 

a M = Ac^ + Bq^ = (x - y + z - u)c^ + (u + z)q jl , 
b^ = + q^ = Ac^ + {B + 1)^, 
h ll =a ll + b ll = 2Ac^ + (2B + 

A = I 2 - [x(k - c) 2 + y(k + c) 2 + z(k + q- c) 2 + u(k + q + c) 2 ] 

= 2(c • q){u 2 — u — xu + yu + z — xz + zy — z 2 } + q 2 {u 2 — u + z 2 — z + 2uz} 

+c 2 {x 2 +y 2 +z 2 +u 2 -x-y-z-u - 2xy + 2xz - 2yz - 2xu + 2yu - 2zu}. (A14) 

Integrating over the loop momenta I, we get 

6e 2 i f 1 

iIii V (q, c) = =— / dxdydz S(x + y + z — 1) 

n 6 Jo 

x j-rT A + c^(^A 2 - H) + v ^c 2 (- 7 -A 2 - 1) + ^(yB 2 + ^-B + I) + r/V(-^ 2 ~ ^ ~ ^)} + f inite 

2 ■ 

= W - <f<f) + finite. (A15) 

It is interesting to observe that the divergent part is not modified by the background chiral gauge field although it 
makes the expression much complicated in the intermediate stage. Contributions from the background chiral gauge 
field turn out to be cancelled exactly in the polarization contribution. It is quite laborious to check this cancellation. 

We summarize leading divergent contributions of fermion one-loop self-energy, one-loop gauge-fermion vertex, and 
gauge-boson one-loop self-energy as follows 

2 ■ 2 ■ 

e % e % 

-iEiQ*,c) = g^2^+ g^2^T5 + finite, (A16) 
T^(p,p', c) = -^7" + finite, (A17) 



As a result, we obtain 



ill^(q,c) = - — Wq 2 - q»q») + finite. (A18) 



e 2 

5$ = ~^T e + finite, (A19) 
e 2 

S c = —^-+ finite, (A20) 
e 2 

S e = —^T e + finite, (A21) 
e 2 

S A = -— r + finite. (A22) 
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3. Renormalization group equations 

We would like to point out that the Ward identity of Z^, = Z e is satisfied, guaranteeing the gauge invariance. We 
emphasize that this serves quite a nontrivial check for our renormalization group analysis, where complex dependencies 
for chiral gauge fields are all cancelled to give rise to the Ward identity. Recalling our perturbative analysis in the 
one-loop level, the satisfaction of the Ward identity implies that the renormalization group equation for the coupling 
constant does not change, compared with the case in the absence of the background chiral gauge field. On the other 
hand, the renormalization group equation for the chiral gauge field is an essential point of our study. 

The beta function for the chiral gauge field is given by 

dc 

/UM)=/^r- (A23) 

Considering that the bare quantity cbv — Z c c v is independent of the scale parameter of /i, we obtain the renormal- 
ization group equation for the chiral gauge field 

d , dM c de d , 

0=- lnc B „ = — S.- + - lnc„, A24 

aln^i de dm^i dm^i 



where 



M c =lnZ c = ± ^tA = -^ + 0(e<) + o{ l_ y (A25) 



e" e f 

n=l 



In the one-loop level we obtain mi(e,c) = Inserting = /3 e (fi) — ee into the above expression, we reach the 

following formula 



- i (-^2 + )\ + ee ) + ^/W (A26) 



Renormalizability guarantees the cancellation in higher negative orders. As a result, we obtain 

dc v e 2 

dp, ' 47T 2 



M = V-^ = --^c„ + 0(e i ). (A27) 



4. Low temperature behaviors for background chiral gauge fields 

Solving the renormalization group equation for the coupling constant 

e 2 M = f° — ■ (A29) 

Substituting this solution into the renormalization group equation for the chiral gauge field 



we obtain 



we obtain 



= In 



4tt 2 ( D ) 

dx 

in d [a^ 1 +lnD] -x 

1 

dx — r 

x — a D 



(A31) 
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where old — j^i is the fine structure constant at the cutoff scale. As a result, we find 



C„(/i) = C V {D) OLD 



1 V 



'D 



(A32) 



Appendix B: Quantum Boltzmann equation approach in the presence of both chiral anomaly and gauge 

interaction 

1. A formal development of the quantum Boltzmann equation in the presence of the topological E ■ B term 

We start from the quantum Boltzmann equation for a steady state [24] 



dG < {p 1 uj) t . ^dG < {p 1 uj) fdZ<(p,uj) d$tG ret (p,uj) d^G ret (p,uj) <9£<(p,w 



. .0U^{p,LJ . r. 
+ rp dco - p -Y 



dp 8lj I dcj dp duj dp J 

= -2T(p, ^)G< (p, w) + £< (p, w)i4(p, w). (Bl) 

G < {p 1 oo) is the lesser Green's function, regarded as a quantum distribution function, where p and uj represent 
momentum and frequency for relative coordinates, respectively. O denotes the derivative with respect to time t. 
S < (p, oj) and G ret (p,cu) indicate the lesser self-energy and the retarded Green's function, respectively, where 3ft is 
their real part. The right hand side introduces collision terms, where F(p,w) and A(p,w) indicate the scattering rate 
and the spectral function. 

r and p are governed by semi-classical equations of motion [7] , given by 



de v 

r = d +PX p ' 



p = eE+-rxB, (B2) 
c 



where Q p represents the Berry curvature of the momentum space. Solving these equations, one obtains 

r = (l + -JB ■ fi p ) {v p + eExSl p + Ul p -v p B}, 

\ + -B-npj 1 {eE+ -v p x B + — (E- B)fi p j. (B3) 



An essential point is the presence of the E ■ B term in the second equation, imposing the Adler-Bell-Jackiw anomaly 
[12]. 

Inserting these equations into the quantum Boltzmann equation and performing some algebra, we obtain the 
following expression 

/ e \ — 1 6 / dG^ \ / 6 \ — ^ /e \ dG^ 

(l + -B-n p ) -v p .(Bx—) + (l + -B.n p ) (- Vp xB).(eExn p) — 

-[A(p, ^)] 2 (-^^) (l + ~ C B ■ «p) ~\eE + € -{E ■ B)n p } ■ {v p + f -(^ p ■ v p )b}t 

= -i[2TG< - S < A], (B4) 

where the argument of (p, ui) is omitted for simplicity. 
The lesser self-energy is given by [25] 

E<(p,w)=V / dv 1P ^- 2 ^D a {q,u){[n{u) + l]G < {p + q,uj + u)+n{u)G < {p + q,LO-u)\. (B5) 
, Jo ™> I J 

Here, we consider gauge interactions for example. Thus, D a (q, v) represents the Green function of gauge fluctuations. 
n(y) is the Bose-Einstein distribution function. One can replace the gauge-boson propagator with some other types 
of fluctuations such as phonons, spin fluctuations, and etc. One may consider the diffusion-mode propagator for weak 
anti-localization, where the form of its vertex should be changed, of course. 
We write down the lesser Green's function in the following way [25] 



G < 



(p,co) = if( u )A(p, U )+i(-^l)A(p,u)v p ■ A(p,w), (B6) 
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which consists of the equilibrium part (the first term) and its correction term (the second term). We call A(p,u>) 
"vertex distribution function" although it sounds somewhat confusing. f(ui) is the Fermi-Dirac distribution function. 

Inserting this ansatz into the quantum Boltzmann equation with the expression of the lesser self-energy and per- 
forming some straightforward algebra, we obtain 



i 



i 1 + -B ■ n, 



-b n 7 



v p x B) ■ (eE x Q p ) 



s 2 
c 

dv 



A(p,w)(l + -B • n p ) 2 [eE +^{E- B)O p } • \v p + ^(n p • v p )B}r(p, cj) 
2T(p,u)v p -A(p,u>)- V / 



pxq 2 



ZD a (q, v){[n(i>) + f(u + v)]A(p + q,uj + v)v p+q ■ A(p + q,u + v) 
-[n(-v) + f(us - v)]A{p + q, u - v)v p+q ■ k{p + q,uj - (B7) 



where we have used the following relation 



21>,u;) = £ 



dv 



pxq 



m 



~SD a {q, v)\[n{v) + f(u + v)]A{p + q,u + v)- [n{-v) + f{uo - v)]A{p + q,u - i/)}. 

(B8) 



Writing down the quantum Boltzmann equation in terms of components and focusing on dynamics near the Fermi 
surface, we reach the following expression for each component, 



AIM 



i(i + -B-n F ) ^A y F (co)-i(i + -B-n F ) 

V c / mc V c / 



mc 



A F ( W ) 



= -i(l + ^B • ft F ) |^B s (e£ x J7 F ) 2 - ^B z (eE x n F ),,} 

— 2 2 2 

-A(p F ,w)r(p F ,w)(i + ^B-n F ) {e£ x + ^(B-B)rj F + ^(i + ^B-n F )(B-B)rj F }, (B9) 



^l-i(l+ e -B. fi F ) _1 — A F (w) + i(l + ■ H F ) "^A F M 
Tt r (w) ' c / mc c / mc 



= + -B ■ fi F ) 2 {- e -B x {eE x fl F ) z + -B z {eE x fl F ) x } 

— 2 2 2 

-A(p F ,uj)T(p F ,u)(l + e -B ■ Sl F ) {eE y + ^(E- B)Q F + ^ (l + -B ■ n F ) [E ■ B)ft F }, 



(BIO) 



and 



(i + -B-n F ) '^H-jfi 

V c / mc V 



c J mc FK ' 



Ttr(w) 

= + ■ n F ) {--B y (eE x fl F ) x + *^B x (eE x fi F )j,} 

— 2 2 2 

-i4(pjr,w)r(p F ,w)(l + -B • n F ) {e£ z + ^-(B • B)fi F + + ^B • H F )(B • B)fl F }, 



(Bll) 



where the transport time is given by 



E 



p F xq 



m 



ZD 



a {q,v)(l - cos 6>){[n(z/) + /(w + v)\A[p F + q,u + v) 



T tr (uj) 

-\n(-v) + f(u - v)]A{p F + q,w - 
We note the 1 — cos 6* factor in this expression, which extracts out back scattering contributions. 



(B12) 
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2. Current formulation 

It is natural to define a current in the following way [12] 

J=-e 



7?E / ■S^{l + lB-^ p y 1 {v p + eEy.n p + - c {n p -v p )B}[- l G<(pM\- (B13) 

iui 



We note the r term in the integral expression. 

Inserting the ansatz for the lesser Green's function into the above expression, we obtain 



)f(u)A(p,u) 



-4 ? / 4k t 1 + -c B • n ^ ~^ vp + > p - vp)B } (-^) ^ • A(p ' w) - (B14) 



Then, each component is given by 
2 - / d3 P 



-e 



-e 



(27Tft) S 



-^E / TOT^ + ~ C B ■ ^"^{K^A^a;) + (^) 2 fi£A,(p,u;) + (^) 2 ^A,(p,c)}(-^)yl(p, w ), 



(B15) 

and 

J . = 4E / 7^(l + ^^ P )"(^-M P )/HA(p,c,) 

+4E / 7^3 (! + I 13 ■ ^pY'-B^vlfniKip^) + (vtfnvAyfauj) + (v;fn p A z (p,u)}(-^±)A(p,u). 

ioj ^ ' 

(B16) 



3. Longitudinal magneto-transport 

Solving the quantum Boltzmann equation in the unconventional setup of B = B x x and E = E x x, we find 

Af.(w) = -eA{p F ,u)^PrE x , (B17) 



a^H = rue (-n* + n^ I f^L 

Ttr(w) 2 2 £C / \ 

-A( PF ,u)- + e -(l + ^)}(n* + n . ^Z^L )E,B„ (B18) 
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and 

A fM = me- (l* + ^ T f^K 

-^.«) 7 W + 7(1 + ^,n%)}(n- F - n v ^Z^L ) £ ,B„ ( Bi9) 

where the "cyclotron" frequency associated with the B x field. We notice that there arc E B = E X B X 

terms, which are topological in their origin. 

Inserting these vertex distribution functions into the current formula, we obtain a rather complicated expression 
for the x— component of the current, 

e 3 1 ^ f d 3 p 1 2 / 2 r fr H 

+-^2.7 M 3 l + §W p) H^-^— ^B-^. 

me 3 1 /" 1 / y \2 / df(u>)\ ., . 

-— ^E/ (2 , fi )3 1 + ^ F «) 2 (-^r)^^) 

B x n*^ 2/ , a/H^ r/1 ^ ^ i2 t^)((°f) +^F fi FTTfB^J F 



me 3 1 v-^ /" <i 3 p 1 , ..,/ df(u>)\ ., . 

— ^E y ( J)»i + f B, n ^ ^ (-^r) A ^") 



+^E/ (^irp^^) (-^r)[^.«)] (l + f ^)V torM]2 (B20) 

Expanding the above expression up to the second order for the Berry curvature and keeping only even-power 
contributions [12], we obtain 

me 3 1 v-^ f d 3 p . v .~( df(u>)\ A . , (Q}L) 2 uj*T tr (uj) 



- t^WW + ^W ^^. lg* - 2C^^ 7!f^ff 2 B, E ,. (B21) 

The first term is also the conventional contribution near the Fermi surface, but there is no dependence for magnetic 
fields. This is certainly expected because the magnetic field is in the same direction as the electric field. On the other 
hand, the second contribution originates from the topological E ■ B term. The third term is also anomalous, which 
results from the Berry curvature but not from the E ■ B term. 
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4. Discussion 

The longitudinal magnetoconductivity is 

„(«., D . CN F eV r r tA T) + g^ ^^ g - ^^ Trfg^ ^ (B22) 

If we limit our discussion on low magnetic fields, we are allowed to neglect the last contribution. Then, the above 
expression can be rewritten as follows 

a L (B x ,T) = (l + C w B 2 x )a n (T), (B23) 

where <r n (T) — CNpe 2 VpT t r(T) is the normal conductivity and Cw — 2(C'/C)(e 2 /c 2 ) is a positive constant. 

Our proposal is to replace the B 2 , term with [c(T)] 2 , where c(T) represents the distance between two Weyl points. 
Then, the final expression for the "longitudinal" conductivity becomes 

a L (T) — ► (1 + /C[c(T)] 2 K(T), (B24) 

where K. is a positive numerical constant and the normal conductivity is determined by intra-scattering events at one 
Weyl point. 
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